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For a two-state quantum object interacting with a slow mesoscopic interacting 
spin bath, we show that a many-body solution of the bath dynamics conditioned 
on the quantum-object state leads to an efficient control scheme to recover the lost 
quantum-object coherence through disentanglement. We demonstrate the theory 
with the realistic problem of one electron spin in a bath of many interacting nuclear 
spins in a semiconductor quantum dot. The spin language is easily generalized to 
a quantum object in contact with a bath of interacting multi-level quantum units 
with the caveat that it is mesoscopic and its dynamics is slow compared with the 
quantum object. 

Coherent superposition of states of a quantum object is the wellspring of quantum prop- 
erties and key to quantum technology. Decoherence of a quantum object results from the 
entanglement with an environment by coupled dynamics |liy,|3[. Amelioration of decoher- 
ence becomes important in anysustained quantum process. Different types of amelioration 
include dynamical decoupling |^, la, |6[ , decoherence- free subspace |7| , quantum error correc- 
tion (for a review, see 8|), and feedback control [9|. 

We offer an alternate approach to the restoration of coherence based on the theory that 
control of the quantum object can direct the quantum evolution of the bath to disentan- 
gle the object from the bath. The operation resembles the spin echo schemes |lO| but it 
removes the pure decoherence due to bath interaction dynamics as well as the inhomoge- 
neous broadening effect. The key is that the environment is effectively a mesoscopic system, 
i.e., the number of particles N is small enough for the timescale of the quantum object 
decoherence to be much smaller than its energy relaxation time Ti while large enough for 
ergodicity, specifically for the Poincare period to be effectively infinite as compared to Ti. 



*Presently at Department of Physics, The University of Texas, Austin, Texas 78712-0264 
^Presently at Department of Physics, The Chinese University of Hong Kong, Hong Kong, China 



The theoretical demonstration of coherence restoration uses one electron spin in a semi- 

digmatic 
III and 



conductor quantum dot of many (A^ ~ 10 ) nuclear spins, which serves as a paradigmatic 
system of a two-level system in a bath of interacting spins for decoherence physics 
for spin-based quantum technology 12,ll3|- Electron spin decoherence due to the hyperfine 
interaction with the nuclear spins has been much studied (see jlJ| and references therein). 
Theories of the effect of interaction between nuclear spins on the electron decoherence have 
recently appeared J3|, 15| . Our theory of coherence recovery by disentanglement is based on 
the previous finding 3] that the mesoscopic bath of slow dynamics is well described by a 
simple pseudospin model for the particle pair interaction in the bath. 

We start by defining the decoherence process in a simple quantum prescription. The 
model for the coupled spin-bath system is a localized electron of spin 1/2 coupled to a bath 
of finite A^ mutually interacting nuclei with spin j in a magnetic field. The isolation of the 
electron spin plus the meso-bath, in the timescale of the total processing duration Tp, from 
the rest of universe arises out of their weak coupling with the outside. The initial state of the 
electron spin, \{p'^{0)) = C+|+) +C_|— ), is prepared as a coherent superposition of the spin 
up and down states |±) in an external magnetic field. The state of the total system of the 
spin plus bath at that instant forms an unentangled state, |^I^(0)) = \ip^{0)) ®\J'). It evolves 
over time t to the entangled |^(t)) = C+(t)|+) 0\J^{t)) + C_(t)|-) ® \J'{t)) where the 
bath states \J''^(t)) are different. The electron spin state is now given by the reduced density 
matrix by tracing over the bath states, p^_^/(t) = C;,{t)C^{t){J''' {t)\J''{t)), a,a' = ±. The 
environment-driven shifting between p^ _,_ and pi _ is longitudinal relaxation. Either off- 
diagonal element gives a measure of the coherence of the spin. The longitudinal relaxation 
contributes to the decoherence. When this contribution is removed, the remaining deco- 
herence is called pure dephasing. For applications in quantum technology, the longitudinal 
relaxation time Ti can be made much longer than the processing duration Tp by a choice 
of system and of the electron Zeeman splitting much larger than the dominant excitation 
energies in the bath and the spin-bath coupling strength |l2l, UM, ll2|- When the cause of 
spin flip is removed, the reduced Hamiltonian of the whole system is in the form diagonal 
in the electron spin basis, H = |+)(+| ® H^ + |— )(— | Cg) H^. The bath evolves under the 
Hamiltonians H^ into separate states \J'^{t)) = e^"^^ '^\J') depending on the electron basis 
states |±). Pure dephasing is then measured by £^_(t) = \{J'\e^^ t^-iH ^^j'^y The elec- 
tron spin coherence may be restored by exploiting the dependence of the bath dynamics on 



the electron spin states to make the bifurcated bath pathways intersect at a later time, i.e., 
I<^^(^)) = l<^ (^))5 leading to disentanglement. 

At temperature (~ lOmK— IK) ^ the nuclear Zeeman energy ujn (~ mK) ^ nuclear spin 
interaction (~ nK), the nuclear bath initially has no off-diagonal coherence and is described 
by Ylij Pj\'^){J'\ where \J) = 0„ |j„), in is the quantum number for J^, the component of 
the nth nuclear spin along z (the magnetic field direction), and Pj gives thermal distribution. 
The essence of electron decoherence is contained in the consideration of each pure bath state 
\J) and later the ensemble average over J is included. The transverse interaction Jn-Jm 
between two bath spins creates the pair-flip excitation, |jn)|jm) -^ |jn + l)|jm — !)• We sort 
out all such elementary excitations from the "vacuum" state \J) and denote each by the 
flip-process of a pseudospin 1/2 indexed by k: \ |fc) ^ \ Ik), characterized by the energy cost 
±ii^fc + -Dfc and the transition matrix element ±Ak + B^ depending on the electron |±) state. 
±-Efc (~ -^/^) is from the longitudinal interaction of the form S^J^ between the electron 
spin S^ and each of the two nuclear spins [2i|. Ak (~ A^/N'^Vl, Vl being the electron Zeeman 
energy) is the extrinsic coupling between two nuclear spins mediated by virtual spin-flips of 
the single electron Isf. Its dependence on the number of particles in the bath signifies its 
mesoscopic nature. B^ is due to the transverse part of the intrinsic (e.g. dipolar) interaction 
between the nuclear spins with strength ~ b. The extrinsic interaction A^ couples any two 
spins in the meso-bath, as opposed to the finite-range intrinsic interaction B^.. D^ (~ b) is 
due to the longitudinal part of the intrinsic nuclear interaction. 

In the nuclear bath with the descending order of parameters, Q ^ Un ^ A/N ^ b, the 
bath dynamics is slow and the density of pair-flip excitations created from the "vacuum" 
state \J') is much less than unity in timescale of interest J3|. The excitations are almost 



always spatially separated, leading to the pair-correlation approximation 3|. 



Iq which treats 



pair-flips as independent of each other. The bath is then driven by the effective Hamiltonian 
derived from the first-principles interactions pf, 

^^ = E^^^-Eh^^J2, (1) 



where ct^ is the Pauli matrix for pseudospin k driven by a pseudo-magnetic field h^ = 
{2Bk ± 2Ak, 0, Dk ± -Efc) depending on the electron |±) state. 

From the justification that correlations of more than two spins are negligible [3|, we derive 



the restrictions which the decoherence timescale places on the size of the bath N, given by 
N%'^A'^ < 1 < niin (y/N, N^^Q^A'^] , to be estabhshed below. The upper bound for 
A^ distinguishes the bath from a macroscopic system. It comes from the dominance of the 
pair correlation in the interaction dynamics of the bath spins over the correlations of more 
than two particles due to the intrinsic interaction. The lower bound, N'^b'^Q^A^^ ^ 1, is 
by a similar consideration but due to the extrinsic interaction of the bath spins. The lower 
bound y/N ^ 1 simply signifies the necessary statistics for decoherence. In the case of the 
electron spin in a GaAs quantum dot, the theory is well justified for 10* > N > 10^ which 
covers quantum dots of all practical sizes. 

The theory of the interacting nuclear spin dynamics dominated by the pair excitation in 
the form of the pseudospin evolution leads to a simple physical picture of coherence decay 
and restoration. The initial unpolarized bath state \J') = (^„ |j„) can be replaced by the 
pseudospin product state <S)k\ Tfe)- Each pseudospin, representing a nuclear spin states 
pair, initially points along the pseudospin +z-axis and then precesses about the pseudo- 
magnetic field h^, \^p^) = e^^^'^'"'''^] tfc), depending on the electron |±) state. Thus, the 
electron spin coherence is measured by the divergence of the pseudospin paths, £^_(t) = 
rifc \{i^k\^k)\ — 6"^'='''=/^. 6k = a/1 — \{'ipk\'ipk)\'^ is the geometric distance between the two 
conjugate pseudospin paths on Bloch sphere. 

Now we examine the consequences of the pseudospin echo. A fast vr-pulse applied at 



t = r to flip the electron spin 



2l| would cause the pseudo-spin evolution 



|^±(t)) = e-|'^^*'=(*-^)e-thfe -fe-| I). (2) 

To find out how to control the decoherence, we neglect for the time being the diagonal nuclear 
spin interaction Dk, which contributes to the same component of the pseudo- magnetic field 
as Ek but much smaller. Because the pseudo-fields dominated by the extrinsic nuclear 
spin interaction, h^ = ±{2Ak,0, Ek), invert exactly into each other, disentanglement of the 
electron spin from the affected bath spin pairs follows at 2r as in the classic spin echo to 
remove the inhomogeneous broadening effect. The pseudo- fields dominated by the intrinsic 
interaction, h^ = {2Bk, 0, ±Ek), do not exactly invert under the influence of the electron spin 
flip and the resultant pseudospin trajectories are illustrated in Fig. ^a). Disentanglement 
from the affected nuclear pairs and, hence, recovery of the electron spin coherence occurs. 




"~^\ r /tX 


'; (b) 


LS! 



10 20 30 40 

f(ns) 


lES^^TXSTliXI 



^0 10 20 30 40 50 60 70 80 

FIG. 1: (a) Precession of the conjugated pseudo-spin vectors under single-pulse control. Red (blue) 
trajectory denotes when the electron state is |-|-) (|— ))• (b) Electron spin coherence £^ _{t) under 
the control of a single flip pulse applied at r = 17 ^us (indicated with the blue arrow). The revived 
coherence peaks at v2r (indicated by the green arrow) while coherence at the conventional spin 
echo time 2r (indicated with the black arrow) is negligible, (c) Contour plot of the electron spin 
coherence C^ _ under single pulse control as a function of time t and the pulse delay time r 
(indicated by the left tilted line) . The restoration of the coherence is pronounced at \/2t whereas 
no coherence peak is visible at the conventional echo time 2r (indicated by the right tilted line). 
The horizontal line is the cut for the curve in (c). (d) Electron spin coherence with a sequence of 
TT-rotations (indicated by blue vertical lines) at intervals of r = 10 /US. The numerical evaluation 
is performed on a GaAs dot with thickness d = 8.5 nm in growth direction [001] and lateral Fock- 
Darwin radius ro = 25 nm, on the large- A^ side of the mesoscopic regime where the intrinsic nuclear 
interaction dominates |3|. Bext = 10 T along the [110] direction. The electron ^-factor is —0.44. 
The random choice of the initial bath state {J) = (^^ |j„) from a thermal ensemble at temperature 
T = 1 K makes negligible difference to the results. 



by the rotation kinematics, at time V2t, distinct from the classic echo. Fig. ^b,c) give the 
computed results for the electron spin coherence in a dot of 10^ nuclear spins in GaAs which 
reveals the coherence recovery after a flip of the electron spin at a range of values for r, even 
after the coherence has apparently vanished. 

Furthermore, the coherence may be restored by a sequence of electron spin flips. For 
example, with a sequence of vr-pulses evenly spaced with interval r, the disentanglement 



from the bath will occur at \Jn{n + 1) r between the nth and the (n -|- l)th pulses, as 
illustrated in Fig. ^d). Consider the correction from the small term Di^ in the pseudo- 
fields, the residue decoherence, at the disentanglement point ^Jnin -|- 1) r, is measured by 
b\ ~ (EkBkDkT^)'^. Compared to the free-induction decay [2| where 61{t) ~ E^BIt'^, the 
decoherence is reduced by a factor of ~ -D^r^ (~ 10^'* for r ~ 10/is). 

Ensemble average over the mixed bath states is necessary in two scenarios, namely, ob- 
servation of decoherence of an ensemble of quantum objects and observation of a single 




FIG. 2: (a) Pseudo-spin trajectories (projected on x-y plane) driven by intrinsic nuclear interaction 
with two vr-rotations of the electron spin at r and 3r. (b) Evolution of the electron spin coherence 
under two-pulse control with r = 5 //s. The dashed blue line denotes the pure state dynamics part 
C^ _ and the solid red line includes the inhomogeneous broadening factor £^ _ . The blue arrows 
indicate the times of the electron spin flip, (c) Contour plot of the ensemble-averaged electron spin 
coherence under the two-pulse control. The tilted lines indicate the pulse times and the horizontal 
line is the cut for the curve in (b). In (b) and (c), we have artificially set the ensemble dephasing 
time Tg in C^ _ to 0.5 /is, about 100 times greater than its realistic value, to make the echo visible 
in the plot, (d) Concatenated pulse sequences. The (l + l)th order sequence is constructed by 
two subsequent Ith order sequences, with a pulse inserted if I is even, (e) Dependence of echo 
magnitude on the echo delay time ti under the control of the concatenated pulse sequences in 
ensemble measurement. The quantum dot is the same as in Fig. ^ except smaller, d = 2.8 nm and 
ro = 15 nm. The nuclear bath is assumed initially in thermal equilibrium at T = 1 K. 
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24| . The coherence of the elec- 



quantum object repeated in a time sequence |2ll 

tron spin is now p+,-(t) = C*_C+C%^_{t) x Cl_{t), where Cl.{t) = Y.jPje~"^^^^^ is the 

inhomogeneous broadening factor due to the probability distribution Pj of the initial bath 

state \J') (different nuclear bath state may result in different Overhauser energy-splitting 

n 

£j of the electron) [3]. (pjit) = Sj [ti — {t2 — ti) + ■ ■ ■ + (— l)"(t — r„)] under the control 
of a sequence of vr-rotations on electron spin at Ti,T2, . . . , and r^. The coherence factor 
£^_(t) is insensitive, up to a factor of 1/yN <^ 1, to the selection of initial bath state 
I JT") = (^„ \jn) (verified by numerical evaluations), and is taken out of the summation [3|. 

Both the inhomogeneous broadening and the pure decoherence due to the extrinsic bath 
spin interaction are shown to be removed at the classic spin echo time 2r in contrast to 
the unusual recovery time of \^t in the case of intrinsic bath spin interaction. We need 
a pulse sequence to produce a time where the decoherence from all three sources can be 
removed. A solution is a two-pulse control. Fig. Efa) shows that, after a second electron 



spin flip at 3r, the paths of the two pseudospin states corresponding to the electron |±) 
states driven by the intrinsic nuclear spin pair interaction cross again at At, coinciding with 
the secondary spin echo time for the other two causes. This two pulse sequence is well- 
known as Carr-Purcell sequence in NMR spectroscopies jl(j]. The residual decoherence at 
t = 4r is 5^ ~ 16 {EkBk — AkDk) D\t^ . The restoration by two-pulse control of coherence 
in the presence of both pure and ensemble decoherence is demonstrated by the results of 
numerical evaluation shown in Fig. I2fb,c). The electron spin coherence is restored at 4r by 
the second pulse even when the first spin-echo at 2r has completely vanished, illustrating 

n 

the remarkable observation 12^] that the absence of spin echo does not mean irreversible loss 
of coherence. 

The power of concatenation of pulse sequences has been shown in the context of dynami- 
cal decoupling in quantum computation |j|. Similarly, the control of disentanglement of the 
bath states from the quantum object may be enhanced by concatenation. The pseudo-spin 
evolution with the two-pulse control of the quantum object can be constructed recursively 
from the free-induction evolution Uq = g-^hj, ■<TfcT/2^ ^^^ ^j^g concatenation, f/^ = Ui^iUi_i, 
I = 1,2. The process can be extended by iteration. Fig. El^d), to any level, f/^ = e^*^' ■'^fe/2^ 
where 6i is the rotation vector along the axis of rotation through an angle 9i . Disentan- 
glement occurs at ti = 2'r coinciding with the classic spin echo. For small 9f, the recursion 
relation is Of^-^ = 6^ + 6^ ^6^ x 6^. At each iteration, the rotation vectors of the conjugate 
pseudo-spin states have their mean (^0^ + Oj^) /2 increased by a factor of 2 and their differ- 
ence [6^ — 6^) reduced by a factor of Of" ~ 2'6r (deduced by induction from 9f ~ 26r). 
The decoherence is reduced by an order of b^r^ at ti for each additional level of concatena- 
tion till saturation at the level /q ~ ~ log2(^'^)- Hence, the coherence echo magnitude scales 
with the echo delay time according to exp (— (rj/Ti)^'"*"^) as shown in Fig. Efe). Our result 
shows the protection of electron spin coherence by pulse sequences with interpulse interval 
up to ~ 10/is 

In conclusion, we note that our scheme of restoring the coherence depends on the pure 
decoherence being driven by the interaction in the spin bath and by the domination of the 
bath pair excitation in the slow bath dynamics. The pulse sequence design is borrowed from 



the dynamical decoupling schemes in NMR spectroscopies 10| and in quantum computation 



a 



but the disentanglement method aims directly at the bath dynamics. Our method seeks 
not to eliminate the object-bath interaction by dynamical averaging, but to disentangle 
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by controlling the quantum object to maneuver the bath evolution. Thus, elimination of 
coupling between the quantum object and the bath is not a necessary condition for their 
disentanglement, as illustrated by coherence recovery at t = \Jn{n + l)r where effective 
object-bath interaction does not vanish even in the first order of hyperfine coupling A/N . 
Direct observation of coherence echoes at such magic times is possible with the narrowing of 



inhomogeneous distribution by measurement projection [26[. The control of bath spins may 
well develop into a valuable addition to the collection of armaments of coherence preservation 
for quantum information processing. 

We acknowledge support from NSF DMR-0403465, ARO/LPS, and DARPA/AFOSR. 



[I] W. H. Zurek, Rev. Mod. Phys. 75, 715 (2003). 

[2] E. Joos, H. D. Zeh, C. Kiefer, D. Giulini, J. Kupsch, and I.-O. Stamatescu, Decoherence and 

the Appearance of a Classical World in Quantum Theory (Springer, New York, 2003), 2nd ed. 
[3] M. Schlosshauer, Rev. Mod. Phys. 76, 1267 (2004). 
[4] K. Khodjasteh and D. A. Lidar, Phys. Rev. Lett. 95, 180501 (2005). 
[5] L. Viola and S. Lloyd, Phys. Rev. A 58, 2733 (1998). 
[6] L. Viola and E. Knill, Phys. Rev. Lett. 94, 60502 (2005). 
[7] L. M. Duan and G. C. Quo, Phys. Rev. Lett. 79, 1953 (1997); P. Zanardi and M. Rasetti, 

Phys. Rev. Lett. 79, 3306 (1997); D. A. Lidar, I. L. Chuang, and K. B. Whaley Phys. Rev. 

Lett. 81, 2594 (1998). 
[8] M. A. Nielsen and I. L. Chuang, Quantum Computation and Quantum Information (Cambridge 

University Press, Cambridge, 2000). 
[9] F. Buscemi, G. Chiribella, and G. M. D'Ariano, Phys. Rev. Lett. 95, 090501 (2005). 
[10] C. P. Slichter, Principles of Magnetic Resonance (Springer- Verlag, New York, 1992), 3rd ed. 

[II] N. V. Prokof'ev and P. C. E. Stamp, Rep. Prog. Phys. 63, 669 (2000). 
[12] D. Loss and D. P. DiVincenzo, Phys. Rev. A 57, 120 (1998). 

[13] A. Imamoglu, D. D. Awschalom, G. Burkard, D. P. DiVincenzo, D. Loss, M. Sherwin, and 

A. Small, Phys. Rev. Lett. 83, 4204 (1999). 
[14] W. A. Coish and D. Loss, Phys. Rev. B 70, 195340 (2004). 
[15] W. M. Witzel, R. de Sousa, and S. Das Sarma, Phys. Rev. B 72, 161306(R) (2005). 



9 

[3] W. Yao, R. B. Liu, and L. J. Sham, cond-mat/0508441 (2005). 

[17] T. Fujisawa, D. G. Austing, Y. Tokura, Y. Hirayama, and S. Tarucha, Nature 419, 278 (2002). 

[18] J. M. Elzerman, R. Hanson, L. H. Willems van Beveren, B. Witkamp, L. M. K. Vandersypen, 

and L. P. Kouwenhoven, Nature 430, 431 (2004). 
[19] M. Kroutvar, Y. Ducommun, D. Heiss, M. Bichler, D. Schuh, D. Abstreiter, and J. J. Finley, 

Nature 432, 81 (2004). 
[20] The orbital wavefunction of the electron covers all N nuclei, so the contact hyperfine coupling 

with each nuclear spin is A/N, inverse proportional to the bath size N. 
[21] A. M. Tyryshkin, S. A. Lyon, A. V. Astashkin, and A. M. Raitsimring, Phys. Rev. B 68, 

193207 (2003). 
[22] J. R. Petta, A. C. Johnson, J. M. Taylor, E. A. Laird, A. Yacoby, M. D. Lukin, C. M. Marcus, 

M. P. Hanson, and A. C. Gossard, Science 309, 2180 (2005). 
[23] F. H. L. Koppens, C. Buizert, K. J. Tielrooij, I. T. Vink, K. C. Nowack, T. Meunier, L. P. 

Kouwenhoven, and L. M. K. Vandersypen, Nature 442, 766 (2006). 
[24] A. S. Bracker, E. A. Stinaff, D. Gammon, M. E. Ware, J. G. Tischler, A. Shabaev, A. L. Efros, 

D. Park, D. Gershoni, V. L. Korenev, and I. A. Merkulov, Phys. Rev. Lett. 94, 047402 (2005). 
[25] W. K. Rhim, A. Pines, and J. S. Waugh, Phys. Rev. Lett. 25, 218 (1970). 
[26] G. Giedke, J. M. Taylor, D. D'Alessandro, M. D. Lukin, and A. Imamoglu, quant-ph/0508144 

(2005); D. Klauser, W. A. Coish, and D. Loss, Phys. Rev. B 73, 205302 (2006); D. Stepanenko, 

G. Burkard, G. Giedke, and A. Imamoglu, Phys. Rev. Lett. 96, 136401 (2006). 



10 

Supplementary Information for "Restoring Coherence Lost to a 
Slow Interacting Mesoscopic Bath" 

Pseudo-spin model 

We give the salient points of the method of solving the electron nuclear spin dynamics 
in a GaAs quantum dot. The system consists of an electron with spin vector Sc and A^ 
nuclear spins, J„, with Zeeman energies Q and Un under a magnetic field -Bcxt, respectively, 
where n denotes both positions and isotope types (^^As, ®^Ga and ^^Ga). The interaction 
can be separated as "diagonal" (or "longitudinal" ) terms which involve only the spin vector 
components along the field (z) direction and "off-diagonal" (or "transverse") terms which 
involve spin flips. Because Q is much larger than a;„ and the strength of the direct nuclear- 
nuclear interactions (e.g., nuclear dipolar interaction), the off-diagonal part of the electron 
nuclear hyperfine interaction can be eliminated by a standard canonical transformation, with 
the second-order correction left as the hyperflne-mediated nuclear-nuclear interaction. For 
the same reason, the off-diagonal part of the nuclear-nuclear interaction includes only terms 
which conserve the Zeeman energies, thus excluding the hetero-nuclear terms. The total 
effective Hamiltonian can be written as H = He + H^ + ^^ |±)i7^(±|, with H^ = ^S^, 
H^ = UnJ^ H^ = ±Ha + Hb + Hd± He, and 

rV _ Sr^'(hiOm j+ j- _ \^' A j+ J- (n \ 

n-A — / ^ .^ J'n •-'m — / , ^n,m-Jn "^mi W^J 

Hb = 2_^ Bn,rnJn-Jm (3b) 

Hd = 2_^Dn,rnJnJm (3c) 

He = J2 («"/2) Jn^H ^-^n, (3d) 

n n 

where |±) are the eigenstates of S*^, the summation with a prime runs over only the homo- 
nuclear pairs, the subscript A denotes the hyperflne mediated nuclear-nuclear interaction, 
B the off-diagonal part of the direct nuclear-nuclear interaction, D the diagonal part of 
the direct nuclear-nuclear interaction, and E the diagonal part of the contact electron- 
nuclear hyperflne interaction. The hyperflne energy, determined by the electron orbital 
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wavefunction, has a typ.cal energy scale £„-.„-#- 10» s- fo. a dot with about 
10^ nuclei |l[, where A is the hyperfine constant depending only on the element type. The 
direct nuclear-nuclear interaction, which is "short-ranged" (referred here as decaying no 
slower than dipolar), has the near-neighbor coupling Bn.m ~ Dn^m ~ 6 ~ 10^ s^^. The 
hyperfine mediated interaction, which couples any two nuclear spins that are in contact 
with the electron and is associated with opposite signs for opposite electron spin states, has 
an energy scale dependent on the field strength, An,m ~ Wn^^^^ ^~^ ^°^ ^^^^ ~ 40-1 T. This 
hyperfine mediate interaction is differentiated from the "short-range" direct nuclear-nuclear 
interaction by the qualifier "infinite-range" . We work in the interaction picture defined by 
He and Hn in which the dynamics are determined by H^. 

The basis set of the bath are eigenstates of Hn'- (S^n bn)- ^^ Eqn. (jHl), Hd and He are 
diagonal in this basis. The off-diagonal terms Ha, Hb are weak perturbations that will 
excite the bath initially on an arbitrary configuration \J) = \ji) ■ ■ ■ \Jn)- The elementary 
excitations in the nuclear spin bath are pair-fiip excitations created by operators J^J~ in 
the reduced Hamiltonian. Starting from any initial nuclear configuration, the evolution of 
the nuclear spin states by these elementary excitations is of the hierarchy as shown in the 
left side of Fig. IHl We can regard the zeroth layer of this hierarchy, the initial state, as the 
'vacuum' of the pair-fiip excitations and layer n corresponds to n pair-fiip excitations have 
been created. The state at time t is a linear superposition of all possibilities: 

Uit)) = Cjit)\J) + Y,C^At)JlJn\J) + E Ci^P,n^At)JtJvJlJn\J) + " " " • (4) 

m,n l,p,m,n 



Exact hierarchy 



Model hierarchy 




JU~J:K\j)* 



jxJ:j;\3) 



j:j;j:j-j;j;\j) 



2 pair-flips 3 pair-flips 



(7;(T;|a)' 




1 pseudo 2 pseudo 3 pseudo 



■XI a) ^1^ : 



^l^l^U^) 



FIG. 3: Hierarchy of the nuclear spin dynamics. 
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where the summation over the indexes m,n,l,p, . . . are defined such that \J'), J^JnlJ')) 
J^ J^ J^J^\J) , . . . denote different eigenstates of li^^ orthogonal to each other. 

We solve this dynamics in the nuclear spin bath based on a pseudo-spin model as described 
below. We have, on sites labelled ra = 1, . . . , A^, the nuclear spin states |j„) with — j < j„ < j 
for nuclei of spin j. As the elementary excitations are pair dynamics driven by J^Jn i ^^ 
first sort out the pair states \in)\im)- These pair states are divided into three categories: 

1. Down States: A down state \jn)\jm) has a partner |j„ + l)\jm — 1) created by 

JnJm\Jn)\jm) = V(j + jn + 1) {j " jn) V U " Jm + 1) (j + Jm) \jn + l)\jm " 1) 

A down state must have (j„ < j,jm > —j)- There are (2j)^ down states for each bond. 

2. Up States: An up state \jn)\jm) has a partner |j„-l)|jm + l) created by i^im|j„)|jm)- 
The up state must have (j„ > —j,jm < j)- There are (2j)^ up states for each bond. 

3. Bachelor States: A single pair state has no partners connected by JnJm\Jn)\jm) or its 
Hermitian conjugate, i.e., jn = jm = j or j„ = j^ = -j. 

The Bachelor states may be mapped to pseudo-spin states. Since they are scalar states, 
their Hamiltonian terms will commute with every other operators, they can only contribute 
to the phase factor in the electron spin coherence through the Overhauser field, causing 
inhomogeneous broadening. 

We shall find by an explicit construction that the up and down states can be paired to 
provide states of (2j)^ two-level systems - the pseudo-spins. These states are divided into: 

1. Monogamy States: Each state belongs to only one pseudo-spin although its partner 
may be a bigamist. These states are edge states in that at least one of the two spin quantum 
numbers (j„ or j^) equal to ±j but they cannot both be equal to j or to —j. Half of them 
(4j - 1 states) are down states, |j„ = -j)\jm > -j) or |j„ < j)\jm = j), see Fig.H^a). The 
other 4j - 1 states, |j„ > -j)\jm = -j) or |j„ = j)\jm < j) are up states, see Fig. H^c). 

2. Bigamy States: Each belongs to two different pseudo-spins. They are interior states: 
~j < jn < j and —j < jm < j- There are (2j — 1)^ of them, see Fig. Efb). 

For a bond between two sites (n, m), all possible pseudo-spins excitations are sort out 
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FIG. 4: Illustration of construction of pseudo-spin states from a pair of nuclei of spin 3/2. The 
solid dots show state \jn)\jm)- The solid arrows between lines show to which state the operator 
Jn Jm would lead and the dotted arrows show to which state JnJm would lead. (a)Monogamy 
states which are mapped to pseudo-spin down: \u,v,[). (b)Bigamy states which are mapped to 
two pseudo-spins with one up and one down respectively: \u,v — 1,|) \u,v,[) (c) Monogamy 
states mapped to pseudo-spins up: \u^v — 1, j). 
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and labelled by using two numbers u{jn,jm) and v{jn,jm) for -j < jn,jm < j, 

U = jm+ in (5) 

^ 7)\Jn Jm) 

The construction follows from the angular momentum addition of the two sites and is illus- 
trated for i = 3/2, 

u= -3 -2 -1 1 2 3 

3 
2 



111 
2 2 2 



v= 



(6) 



1 _i _i 

'2 2 2 



When u = ±2j, the two-spin states are 

\jn)\jm) = \j)\j) or |-J)|-J). (7) 

Both are Bachelor states. For a u value equal to or between — 2j -|- 1 and 2j — 1, the v state 
at the bottom of the ladder is chosen to be the spin down of the doublet and one rung above 
to be its partner. Thus, the bottoms of the ladders yield the down monogamy states and 
the top rungs the up monogamy states. In between, the states are the bigamy states. 

For any initial configuration \J'), the relevant set of pseudo-spins Qj is determined by 
examining every possible nuclear spin pair {m, n). Each pair will contribute 0, 1 or 2 pseudo- 
spins if \jm)\jn) is in the Bachelor state. Monogamy state or Bigamy state configuration 
respectively (see Fig. HJ. The many nuclear spin initial state \J') is then replaced by, 

|J)^|ji)...|j^)^|a)=(g)|fca) (8) 

keGj 

where k labels both the nuclear pair (m, n) and the pseudo-spin type (u,v). u =1 or | 
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FIG. 5: Illustration of multi pair-flip excitations in the nuclear bath. We use hollow arrowheads 
to indicate the first pair-flip and solid arrowheads to indicate the second pair-flip. Dotted arrows 
denote the pseudo-spin flip from 'down' state to 'up' state and solid arrows for the inverse process 
(see text), (a) Independent pair flips; (b-e) Various situations of overlapping pair-flips; (f-h) 
Approximations in the independent pseudo-spin model. 



depending on whether k is mapped from up or down Monogamy or Bigamy state. Different 
initial nuclear configurations will result in different sets of pseudo-spins. For a randomly 



chosen initial configuration |jr), the number of pseudo-spins is given by M 



- 2j ■ 
' 2j+l . 



''ZN 



where N is the total number of nuclear spins and Z the number of nuclei coupled to a par- 
ticular nuclear spin by the nuclear-nuclear interaction. For short ranged direct interaction, 
Z ~ 0(10) and for the infinite-ranged hyperfine mediated interaction Z = N. The factor 
(2^)^ arises as the single state, Monogamy state and Bigamy state are contributing 0, 
1 and 2 pseudo-spins respectively. For convenience, when the set Qj is determined from 
the \J'), we redefine the pseudo-spin up and down states, i.e. \k,a) -^ \k, —a), for those 
pseudo-spins in set Qj so that the initial state |5) in this new definition corresponds to all 
pseudo-spins pointing 'up': (^j^, | 1)^. The redefinition is conditioned on the initial nuclear 
configuration \J'). 

Conditioned on the electron state |±), the pseudo-spins are driven by the effective Hamil- 
tonian of the form, 

4t = E^^^-Eh^^^/2 (9) 

k k 

The effective magnetic field h^ on the pseudo-spins, conditioned on the electron spin state, 
are to be determined by reproducing the matrix elements, {J\J^J~H^J:^J:^\J) — {J\H^\J) 
and {J'\H^J^J^\J'), namely the energy cost and transition matrix element for nuclear pair- 
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flips. 

The pseudo-spin model for characterizing the nuclear spin bath dynamics is to approx- 
imate the exact evolution of Eqn. ^ by the independent evolution of all pseudo-spins in 

\m) = (g) i^kit)) = c,im) + $^c,,(t)a+ia) + J2 Ck,MitK&) + ■■■■ (10) 

k£Qj k\ ki,k2 

This pseudo-spin dynamics can be put into a similar hierarchy as shown in the right part of 
Fig. 121 which we will refer to as model hierarchy in contrast to the exact hierarchy. 

With the mapping established for the state (Eqn. (jH))) and the Hamiltonian (Eqn. Q), 
the first two layers of the exact hierarchy will be reproduced exactly by the model hierarchy, 
i.e., there is a one to one correspondence between J^JnlJ') and (3"^|5) with the energy and 
coupling to the initial state \J') (|5)) exactly reproduced. 

Difference between the model and the exact hierarchies arises when more than one excita- 
tions have been created in the system. In Fig.^l we illustrate with the case when two pair- flip 
excitations have been created. If the two pair-flips do not overlap as shown in Fig. Efa), 
their dynamics are then independent of each other and well described by the pseudo-spin 
model. FigE^b-e) illustrate the various situations that the two pair-flips overlap, by sharing 
one or two nuclei. The flip-flop of the first nuclear pair (/, m) changes the spin configuration 
of both nuclear / and m and if a second flip-flop is to take place on pair (/, m) or (n, m) 
or (/,n), it is no longer described by the dynamics of the original pseudo-spins assigned to 
it. Instead, in the model hierarchy by the independent pseudo-spin model, two successive 
flip-flops on pair {l,m) or two successive flip-flops on pair {l,m) and {m,n) respectively are 
shown in FigEl^f-h). Fig. Efg) can be considered as the approximate form of Fig. El^c) and 
Fig. Efh) as that of Fig. Efd). The model hierarchy contains events like Fig. Eff) which is 
absent in the exact hierarchy and events like Fig. Efe) in the exact hierarchy is not contained 
in the model hierarchy. Therefore, on layer 2, the model hierarchy coincides with the exact 
hierarchy in events described by Fig. Efa) and differ by replacing the events of Fig|Sl^b-e) 
with events of FiglH^f-h). The difference in a general layer can be analyzed in the same way. 

By the pseudo-spin model, we are using Eqn. (fTIHl as the bath state at time t for calculating 
physical properties instead of Eqn. (Q. The difference of the exact and model hierarchies 
is estimated below which serves as an upper bound for error estimation (notice that the 
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physical properties of interest are not necessarily changed by replacing events of Fig|Htb-e) 
with events of Fig|Sl^f-h), therefore, this error-estimation is not necessarily a tight bound). 
If n — 1 pair-flip excitations have already been generated, to create the next excitation, we 
have M pseudo-spin to choose from and ~ 2{n — 1)Z of them overlap with the previous 
excitations. Therefore, the probability of having a new excitation without overlapping with 
the previous excitations is given by: ~ ~ m ■ By induction, the probability of creating 
n non-overlapping pair-flip excitations is then given by, 

M -2Z M -2(n-l)Z 
pn ~lx X ^— '— (11) 



exp 
exp 



2Z 2(n-l)Z 



M M 

n{n-l) 2j + l ^ 

N ^ 2j \ 



The second ~ holds if 2nZ <^ M which is always true in the timescale relevant in our study. 
Comparing the model hierarchy and the exact hierarchy, we find from the above analysis 
that they differ in layer n with the relative amount of 1 — p{n). 

Error estimation 

We perform a self-consistent analysis on validity of the pseudo-spin model. Here the 
relevant timescale plays the crucial role in determine the validity of the model. If only the 
first n layers of the hierarchy (the exact one and the model one) are involved, the error in 
the calculated physical properties is bounded by 1 — p{n) '^ 1 ii n^ <^ N. Therefore, in the 
very short time limit, only the first several layers of the hierarchy can be involved and the 
pseudo-spin model gives an almost exact account of the dynamics. To estimate error upper 
bound for the longer time limit, we will calculate the number of layers involved (denoted as 
n) based on the model hierarchy of the pseudo-spin model. If n^/N obtained is small, we 
conclude that n also faithfully reflects the number of layers involved in the exact hierarchy. 
Therefore, the error estimation based on the pseudo-spin model is faithful and any physical 
properties calculated based on pseudo-spin model is also a good approximation since the 
difference from the exact dynamics is small. Otherwise, the approximation is not good. It 
is established below that the condition n^ -C iV is the origin of an upper bound and a lower 
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bound on the quantum dot size A^ for which our theory is able to deal with. The validity of 
the pair-correlation approach (pseudo-spin model) is indeed in the mesoscopic regime. 

At any given time t, an average excitation number NfUpit) can be defined as follows based 
on the model hierarchy, 

iV^,,p(t) = lx5^|Cfc,(t)|' + 2x 5^|C,,fc,(t)|V3x Y. W|Cfci,fc„fc3 !' + ••• (12) 

fci fci,fe2 ki,k2,k3 

and our analysis J2| shows that the layer-distribution of population in the hierarchy is of a 
normal distribution centered at Njnp, i.e. the population is distributed in layers from layer 
Nfiip ~ \jNfiip to Nfiip + \jNfiip. Therefore, the quantity for characterizing the error upper 
bound is of a very simple form: Perrif) = 1 — exp{—Nji^p(t)/N). 

In the pseudo-spin model, Nfiip{t) defined in Eqn. (fT^ has an equivalent expression which 
is more convenient for evaluation: 

NfUpit) = J2\(i\Ukit)\^)\' (13) 



where f/^ (t) = e~*^fe * is the evolution operator for pseudo-spin k. The contribution can be 
divide into two parts: A^/^jp (t) = N^^^ (t) + Nj^^^ (t). N^i-^ is the number of non-local pair- 
flip excitations and N?i- is the number of local pair-flip excitations that have been created. 
Nfiip (i) and A^^jp (t) have very different behavior and we analyze them separately. 
In free-induction evolution, the number of non-local pair-flip excitations is given by. 



1. \ K. / h 



where h^ = \/E\ + 4/1^ and Ma ~ A^^ is the number of non-local nuclear spin pairs. 
Since the evolution of the non-local pair-correlation is completely reversed by the vr pulses, 
Nfiip (t) is also reversed and N^^^ = at each spin echo time. Therefore, A^/)jp (t) does 
not accumulate in the pulse controlled dynamics and we just need to look at the maximum 
value of A^//jp (t) between echoes. For example, in the control with the equally spaced pulse 
sequence, in order to have the coherence well preserved or restored at spin echo time, the 
delay time between successive pulses is limited by r < T^, where T^ — b~^^'^A~^^'^N^^'^ is 
the Hahn echo decay timescale and ~ lO/xs in GaAs [3|. Therefore, A^/ljp at any time is 
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bounded by Nj^^ {0{Th)) in this scenario. The estimate of the bound on A^ imposed by 
iVj|- for the vahdity of our approach is therefore, 

For GaAs fluctuation dot in a magnetic fleld of 10 Tesla, the above condition is well satisfled 
for A^ > 10^. Estimation for other controlling pulse sequences give a lower bound on A^ 
which is in the same order. 

For the number of local pair-flip excitations, we have a similar expression in the free- 
induction evolution, 



I, \ k / I. 



aNbH 



2^2 



where /if = a/-E| + 4i?| and Mb ~ aN is the number of local nuclear spin pairs, a 
(~ 10 in GaAs) is determined by number of local neighbors and the nuclear spin quantum 
number j. In contrast to the non-local pair dynamics, the local pair dynamics is not reversed 
under the influence of the electron spin flip and N?i-{t) accumulates all through the time. 
Nonetheless, it turns out that Nf^-it) < J^k^l^'^ — oi^b'^'t^ holds for all scenarios of pulse 
controls being discussed. Therefore, the condition {Nfi^^Y /N <^ 1 sets an upper bound on 
A^: Na^ [bt) <^ 1, which depends on the time range t we wish to explore. Alternatively 
speaking, {N^i-p^/N <^ 1 sets an upper bound on the time range t we can explore for some 
flxed A^ using the pseudo-spin model. We illustrate this bound using the following two 
examples. 

1. If we wish to calculate the Hahn echo signal using the pseudo-spin model, we shall have 
{aN}?Tlf/N < 1 where the Hahn echo decay time Th ~ h-^I^A~^''^N^/^ 's]. Therefore, 
the upper bound on A^ is given by N'^a%'^A^'^ ^ 1. For GaAs quantum dot, this condition 
is well satisfled for A^ < 10*^. 

2. For bath of an intermediate size in the allowed region of min[viV, N%'^VL'^A^^] ^ 
1 ^ a^ N'^lP' A""^ , e.g., a quantum dot of typical size A^ ~ 10^ — 10^ in our problem, 
{aNh'^t^Y / N ^ 1 is satisfled for a much longer time range t ~ lOT^ ~ lOOyUS. 

In summary, for the pair-correlation approximation (or pseudo-spin model) to be valid, 
nuclear spin dynamics of local pair-flips imposes an upper bound on A^ while nuclear spin 
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dynamics of non-local pair-flips imposes a lower bound. Within this mesoscopic regime, 
the pair-correlation approximation is well justified. The error estimation is based on char- 
acterizing the difference in the Hilbert space structure of the exact dynamics and that of 
the pseudo-spin model and assuming this difference has a full influence on the electron spin 
coherence calculation. Therefore, the bound is not necessarily tight and it is possible that 
the pulse control methodology developed using the pseudo-spin model have actually a much 
larger validity regime. Investigation is underway. 
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